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Renormalization of the Three-Body System with Short-Range Interactions 
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We discuss renormalization of the non-relativistic three- 
body problem with short-range forces. The problem becomes 
non-perturbative at momenta of the order of the inverse of 
the two-body scattering length, and an infinite number of 
graphs must be summed. This summation leads to a cutoff 
dependence that does not appear in any order in perturbation 
theory. We argue that this cutoff dependence can be absorbed 
in a single three-body counterterm and compute the running 
of the three-body force with the cutoff. We comment on rel- 
evance of this result for the effective field theory program in 
nuclear and molecular physics. 

Systems composed of particles with momenta k much 
smaller than the inverse range 1 / R of their interaction are 
common in nature. This separation of scales can be ex- 
ploited by the method of effective field theory (EFT) that 
provides a systematic expansion in powers of the small 
parameter kR Generically, the two-body scattering 
length 02 is comparable to R, and low-density systems 
with k <C 1/02 can be described to any order in kR by 
a finite number of EFT graphs However, there are 
many interesting systems, such as those made out of nu- 
cleons or of ^He atoms, for which 02 is much larger than 
R. In this case the expansion becomes non-perturbative 
at momenta of the order of 1/02, in the sense that an 
infinite number of diagrams must be resummed. This re- 
summation generates a new expansion in powers of kR 
where the full dependence in ka2 is kept. Consequently, 
the EFT is valid beyond k ^ 1/02, comprising, in partic- 
ular, bound states of size ~ 02. There has been enormous 
progress recently in dealing with this problem in the two- 
body case Q , where the resummation is equivalent to ef- 
fective range theory Q . Ultraviolet (UV) divergences ap- 
pear in graphs with leading-order interactions and their 



resummation contains arbitrarily high powers of the cut- 
off. A crucial point is that this cutoff dependence can be 
absorbed in the coefficients of the leading-order interac- 
tions themselves. All our ignorance about the infiuence of 
short-distance physics on low-energy phenomena is then 
embodied in these few coefficients, and EFT retains its 
predictive power. However, the extension of this pro- 
gram to three-particle systems presents us with a puzzle 
g. Although in some fermionic channels the resummed 
leading two-body interactions lead to unambiguous and 
very successful predictions [^[0], amplitudes in bosonic 
systems and other fermionic channels show sensitivity to 
the UV cutoff, as evidenced in the Thomas |^ and Efimov 
1^ effects. This happens even though each leading-order 
three-body diagram with resummed two-body interac- 
tions is individually UV finite. We will argue below that 
the addition of a one-parameter three-body force coun- 
terterm at leading order is necessary and sufficient to 
eliminate this cutoff dependence. This result extends the 
EFT program to three-particle systems with large two- 
body scattering lengths, including the approach of Ref. 
pOf where pious are treated perturbatively. 

The most general Lagrangian involving a non- 
relativistic boson ip and invariant under small-velocity 
Lorentz, parity, and time-reversal transformations is 
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where the ellipsis denote terms with more derivatives 
and/or fields; those with more fields will not contribute to 
the three-body amplitude, while those with more deriva- 
tives are suppressed at low momentum. It is convenient 
|]lT| to rewrite this theory by introducing a dummy field 
T (called "dibaryon" in analogy to the nuclear case) with 
quantum numbers of two bosons, 

C = ^H^do + + ATtT - -j^iTHi^ + h.c.) 

+ hT''Ti;^iP + ... (2) 
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FIG. 1. (a) Dressing of the dibaryon. (b) Diagrams con- 
tributing to particle/bound-state scattering. 



The arbitrary scale A is included to give the field T 
the usual mass dimension. Observables depend on the 
parameters of Eq. (^) only through the combinations 
Co = gVA = 47ra2/M and Do = -3hg^/A^. The (bare) 
dibaryon propagator is simply a constant i/A, while the 
particle propagator reduces to the usual non-relativistic 
form i/{p'^ — /2M). First, we consider the dressing of 
the dibaryon in Fig. |^(a) at a generic momentum Q. The 
boson loop has a linear UV divergence that is absorbed 
in 3^/A, a finite piece ~ Alg'^Q/AnA^ determined by the 
unitarity cut, and terms suppressed by powers of Q over 
the cutoff A that are subleading and of the same size as 
terms in Eq. (|^) already disregarded. Using the relation 
between g and 02, we see that the finite piece is Qa-ij A. 
and, consequently, has to be resummed to all orders for 
Q^l/a2. As a result, the dressed dibaryon propagator 
to leading order is given by 



(3) 



Attaching four boson lines to this propagator gives the 
two-particle scattering amplitude at leading order, which 
is identical to the effective range expansion at the order of 
the scattering length. Further corrections give the next 
terms in the effective range expansion Q| . 

Let us now consider particle/bound-state scattering. 
All diagrams contributing to this process in leading or- 
der are illustrated in Fig. ||(b). Each of the diagrams 
including only two-body interactions gives a contribution 
of order ^ Mg'^/Q'^ ^ 47rAa2. (The properly normalized 
amplitude is independent of the arbitrary parameter A; 
it appears here only because of our choice of interpolat- 
ing field for the bound state.) The relative size of graphs 
that include a three-body force will be discussed shortly. 

The sum of all the diagrams in Fig. ^(b) can be ac- 
complished by solving the equation represented by the 
second equality in Fig. |l|(b) (|-|7|JlT 



i{p) = K{p, k) + 
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dq K{p, q)- 
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where k (p) is the incoming (outgoing) momentum, 
ME = 3fc^/4 — l/a| is the total energy, a{p = k) is 
the scattering amplitude normalized in such a way that 
as = — a(0) is the particle/bound-state scattering length, 
and 
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The parametric dependence of a{p) on k is kept implicit. 
We are interested in A = 1 for the boson case. This 
equation reduces to the expressions found in Ref. 
when h = 0. Three nucleons in the spin J = 1/2 channel 
obey a pair of integral equations with similar properties 
to this bosonic equation, while the spin J = 3/2 channel 
corresponds to A = —1/2. 

Let us look at the asymptotic behavior of the solution 
of Eq. (^ in the case h = 0. For 1/02 < p < A (but k 
1 /a2) , the integral in Eq. (^) is dominated by momenta in 
the intermediate region 1 /a2 ^ q ^ A and the equation 
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holds up to terms suppressed by powers of l/pa2 and/or 
p/A. The scale invariance of Eq. (^ suggests an ansatz 
of the form a{p) ~ p*, which works if s satisfies 
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The solutions of Eq. (0) come in pairs due to the ad- 
ditional symmetry a(j)) —> a(l/p) of Eq. (|^). For 

A < Ac = ~ 0.4135, Eq. (0) has only real roots. 
However, for A = 1 there are two imaginary solutions 
s = ±iso, with So — 1.0064. Both make the integral in 
Eq. (^) UV finite and are equally acceptable: a{p) is 
given in the intermediate region by a linear combination 
of p^'^° and . Eq. (^ is homogeneous so it clearly 
cannot determine the overall normalization of a, but we 
now see that it cannot uniquely determine the phase ei- 
ther. However, Eq. (^) with finite A has a solution with 
a well determined phase that in the intermediate region 
1/^2 <C p ^ A is, on dimensional grounds. 
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a{p) = A cos(so In — +5), 



(8) 



where 5 is some dimensionless, cutoff-independent num- 
ber. The limit A ^ 00 is not well defined because Eq. 
(^) does not have a unique solution in this limit. (A rig- 
orous proof that Eq. (^) does not have a unique solution 
can be found in p^.) This non-uniqueness comes from 
our idealization of the interactions as point-like. Note 
that subleading contributions from the integration range 
A < g < 00 change phase and amplitude significantly. 
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Numerical solutions of Eq. fl) with fc = for different 
values of A are plotted in FigTg. (Our results agree with 
those in Ref. |0 for the appropriate cutoff values.) We 
observe that indeed the behavior of a{p) in the region 
1 /a2 <C p <C A is given by Eq. and that small differ- 
ences in the asymptotic phase lead to large differences in 
the particle/bound-state scattering length. 
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FIG. 2. a{p) for fc = 0. Full, dashed, and dash-dotted 
curves are for H — and A — 1.0, 2.0, 3.0 x lO^a^^, respec- 
tively. Dotted, short-dash-dotted, and short-dashed curves 
are for A = lO^aJ^ and H — —6.0, —2.5, —1.8, respectively. 

Note that if the series of diagrams in Fig. ^(b) was 
truncated at some finite number of loops one would miss 
the correct asymptotic behavior of a{p) (Eq. (||)) that 
generates the cutoff dependence. This is because Sq (and 
its expansion in powers of A) vanish in a neighborhood 
of A = and the truncation of the series in Fig. [l](b) is 
equivalent to perturbation theory in A. 

This cutoff dependence comes from the behavior of the 
amplitude in the UV region, where the EFT Lagrangian, 
Eq. (|^), is not to be trusted. When the low-energy ex- 
pansion is perturbative, the cutoff-dependent contribu- 
tion from high loop momenta can be expanded in powers 
of the low external momenta and cancelled by local terms 
in the Lagrangian. Thus all uncertainty coming from the 
high momentum behavior of the theory is parametrized 
by a few local counterterms. The present case is compli- 
cated by the non-analytic cutoff dependence of the ampli- 
tude around p = 0. That, however, does not mean that 
the renormalization program is doomed: a three-body 
force term of sufficient strength contributes not only at 
tree level, but also in loops dressed by any number of two- 
particle interactions. This generates non-local contribu- 
tions precisely of the same form as the cutoff-dependent 
terms generated by the two-body force alone. To see how 
that comes about we turn on the three-body force term 
and write h{A) = 2M g"^ H {K) / assuming H{K) ~ 1. 
The asymptotic Eq. (0) becomes 
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where we have set A = 1 for definiteness. For p ^ K the 
term proportional to H becomes important and a(p ~ A) 
has a complicated form. In the range l/a2 <C p <C A, 
however, the three-body force is suppressed by p/A com- 
pared to the logarithm and can be disregarded. Conse- 
quently, Eq. (H) is still correct in the intermediate region. 
The effect of a finite value of H can be at most to change 
the values of the amplitude A and the phase 5, which 
become functions of H. As shown in Fig. ||, this is con- 
firmed by numerical solutions: while different values of 
the three-body force preserve the form of the solution, 
the phase (and amplitude) are changed. If H is chosen 
to be a function of A such as to cancel the explicit A 
dependence, we can make the solution of Eq. (^) cutoff 
independent for all p ^ A. In particular the scattering 
amplitude that is determined by the on-shell value a{k) 
with k ^ I/a2 will be cutoff independent as well. For 
this to be possible A and 5 must depend on the same 
combination of A and H . Thus H{K) must be chosen 
such that 



- solnA + (5(i7(A)) = -sqIhA^, 



(10) 



where A^, is a parameter fixed by experiment or by match- 
ing with a microscopic model. 

We can get a handle on the form of H{h) by consider- 
ing Eq. (^ with two different values of the cutoff A and 
A' > A, whose solutions we denote by a{p) and a'{p). In 
the intermediate region 1 /a2 <C p <C A the equations for 
a(p) and a! {p) will have the same form except for 
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where /i <C A is an arbitrary scale, and we dropped terms 
suppressed by further powers of l/pa2 and p/A. Assum- 
ing a'{p) has the same phase cos(so ln(p/Ai,)) as a{p) even 
for p ^ A', we can make the terms in Eq. ( pT| ) vanish by 
choosing 



^ ^ sin(soln(A/A^) - arctg(l/so)) 
~ sin(so ln(A/A*) + arctg(l/so)) ' 



(12) 



Since K' — K nearly vanishes for all p <C A, a'{p) has also 
the same amplitude as a{p) in the intermediate region. 
That is, with H{A) chosen like Eq. (12), a(p) = a'{p) for 
all values p ^ A (up to terms suppressed by p/A), and 
the on-shell amplitude a{k) for A; <C A will be A indepen- 
dent. Once the parameter A^ is fitted to an experimental 
datum at a certain energy, the energy dependence can be 
predicted. 
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FIG. 3. Three-body force as a function of the cutoff A: 
numerical solution (dots) and Eq.(^2|). 




to be unequivocally within the limits of applicability of 
the EFT; it has a cutoff-independent binding energy of 
B3 = 1.5/Mal 

The value for the ratio azja-i used above is the one 
suggested by the values of 02 = 124.7 A and 03 — 195 A 
given by a phenomenological **He-^He potential |jl^ con- 
sistent with the recent measurement of the dimcr binding 
energy ||l^. Fig. |^ then represents the phase shifts for 
atom/dimer scattering, with an effective range = 71 
A. Similarly, our result for the shallowest bound state 
suggests an excited state of the trimer at B-^ = 1.2 mK. 
Because the integral equations are similar, our arguments 
are relevant for three-fermion systems with internal quan- 
tum numbers as well jl^ . The approach of Ref . ||l^ then 
suggests that our bound-state results would provide a 
reasonable estimate of the triton binding energy. 

In conclusion, we have provided analytical and numeri- 
cal evidence that renormalization of the three-body prob- 
lem with short-range forces requires in general the pres- 
ence of a one-parameter contact three-body force in lead- 
ing order. This opens up the possibility of applying the 
EFT method to a large class of systems of three or more 
particles with short-range forces. 

We thank V. Efimov, H. Miiller, and D. Kaplan for 
helpful discussions. HWH acknowledges hospitality of 
the Nuclear Theory Group and the INT in Seattle. 
Research supported in part by the U.S. DOE (DOE- 
ER-40561 and DE-FG03-97ER41014), the NSERC of 
Canada, and the U.S. NSF (PHY94-20470). 
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FIG. 4. Energy dependence: fccot(5 for different cutoffs 
(A = 42.6, 100.0, 230.0, 959.0 x aj^). 



We also determine H{A) numerically by finding the 
value of H that keeps the scattering length 03 = —a{Q) 
constant for each value of A varying over a large range. 
These values are plotted as a function of A in Fig. || to- 
gether with H{A) given by Eq. (|2|). For illustration we 
used 03 = 1.56a2, but have verified that similar agree- 
ment holds for other values of 03. In Fig. ^ we show the 
corresponding kcotS = ik + a(fc)~^, where 6 is the S- 
wave phase shift for particle/bound-state scattering, for 
several values of A. As argued above, it is insensitive to 
A as long as fc ^ A. The effective range, for example, is 
predicted as = 0.57a2. Note that the three-body force 
discussed here is not the one used in realistic potential 
models where the effective cutoff is at much higher scales. 

These arguments hold for the bound-state problem as 
well because the inhomogeneous terms played no role. 
We have solved the homogeneous equation with the H{A) 
of Fig. 1^. Only the shallowest bound state is large enough 
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